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Abstract 

The theoretical nuclear physics predictions and experimental results on antiprotonic atoms lead 
to different proton and neutron spatial distributions. At CERN SPS energies production of positive 
and negative pions differs for pp, pn, up and nn scattering. These two facts lead to an impact 
parameter dependence of the vr^ to vr^ ratio in ^^^Pb +^08 collisions. A recent experiment at 
CERN seems to confirm qualitatively these predictions. The results for two models considered are 
almost identical. It may open a new possibility for determination of neutron density distribution 
in nuclei. 
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The whole traditional microscopic nuclear physics is based on the assumption that nuclei 
are built of protons and neutrons. The so-called y-scaling analysis seems to confirm that 
this is good assumption 

Electron scattering off nuclei provides direct information about charge distribution (see for 
instance which is closely related to the spatial distribution of protons. The information 
about the neutron spatial distribution is not accessible directly. 

Experimentally one usually determines the difference Ar„p between the neutron and 
proton mean square radii. However, this does not allow one to make a difference between the 
neutron skin or neutron halo model Ij. The radiochemical method applied to antiprotonic 

nnn 

atoms 1^, |5|, |6|| allows one to measure the neutron-to-proton ratio at the peripheries of nuclei. 
This measurement suggests a density profile of the halo type, where half-radii for neutron 
and proton densities are equal, and distributions differ only in their diffuseness. 

Another method of investigating neutron distribution is based on the analysis of X-ray 
spectra of a„tip.oto,nc atomj. It wa. ate =ugge.ted in Q that parity violating electron 
scattering could measure neutron density (for a newer review see p]). 

On the theoretical side the difference between the proton and neutron distributions can be 
obtained in the framework of Hartree-Fock (HE) method (see for example 9]) or Hartree- 
Eock-Bogoliubov method (see for example 10|). Eor heavy nuclei one usually uses zero- 
range nucleon-nucleon interaction of the Skyrme type. Reliable neutron densities can be 
also obtained from a global analysis of intermediate-energy elastic proton-nucleus scattering 

HQ. 

Stable heavy nuclei exhibit an excess of neutrons over protons. A recent experiment 
at CERN for charged pion production in the "^^^Ph +2°^ Ph collision has observed an 
interesting dependence of the ratio / 'i^ '^'^ Eeynman variable xp- Surprisingly, 
the ratio for central and peripheral collisions differs significantly Q]. 

In EigQwe present the ratio Pnij) / Pp{f) as obtained from the HEB densities of "^^^Ph jlfll |. 
The experimental point, taken from Ref.j^l, is a result of interpolation of results obtained for 
other nuclei (the radiochemical method can be used only for a limited set of elements) . The 
grey area shows a direct estimate of Pni^) / Pp{r) obtained from X-ray spectra of antiprotonic 
Pb atoms. This is slightly model dependent and was obtained under assumption of Eermi- 
like density profiles of halo type. At large radii the HEB calculations are consistent with 
these experimental results. 
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In the following we assume that pions are produced in elementary nucleon-nucleon colli 
sions. Without loosing more generality we can write 



{b,XF;W)= N^i3^Mb;W)— ixp;W), (1) 

dxF ^^t;,„ dXF 



where W is energy per binary nucleon-nucleon collision and N^is^n^ are numbers of collisions 
of a given type. Only elementary cross sections for the pp — > vr^ processes are known 



hi 



owing we shall use elementary cross sections -^{xf; W) 
I . We shall leave the problem of the quality of the model 



experimentally. Therefore in the 
calculated in the HIJING model 
for a separate detailed analysis. 

In Fig 121 we show xf dependence of the cross section for different types of collisions. 
We have used = 30 mb to normalize all elementary cross sections. A clear difference 
between different cross sections can be seen. This asymmetry in elementary collisions in 
conjunction with the asymmetry of proton and neutron spatial distributions should lead to 
asymmetries in nuclear collisions which will be discussed in the present note. Only four 
elementary cross sections are independent. The cross sections for the vr" production (right 
panel) can be obtained from those for the vr"*" production (left panel) thanks to isospin 
symmetry. In addition at a;ir > one has (Xpp^^i ~ crp„^^± and cr„„^^± ~ cr„p^^±. 

When calculating the ratio of 7r+ to vr^ cross sections it is sufficient to know the fractions 

naive approach 



faf3 = °fr'' . It is far from obvious how to calculate the fractions fap- In the most 



where Na/i, iV/3/2 are fractions of protons/neutrons in the nucleus Ai and A2, respectively. 

In this simple model there is neither impact parameter nor energy dependence of fap, i.e. 

as a consequence of Eq.((H) the ratio 

,,,,,_ da^^^---\xF,b;W)/dxF 
(xF, b; W) = (3) 

depends only on the Feynman xf and center-of-mass energy. Recent preliminary experi- 
mental results of the NA49 collaboration seem to contradict this simple approach. 

It is often assumed that the dynamics of nuclear collisions is governed by the number of 
binary collisions The number of binary nucleon-nucleon collisions at a given impact 

parameter is proportional to the nucleus-nucleus thickness [3| 

Ta.aM = J dhiT^Asi) TaM - ^) = / dh2TAAs2 - h) T^M) ■ (4) 



In Eq.(jH) we introduced TA^ih) = J dzi pAi(b, Zi), where pa^ is the density function of the nu- 
cleus Ai normahzed to the number of nucleons. Analogously the number of binary collisions 
of a given type in Eq.((H) can be written as 

N^P^Mb; W) = T^aM ■ ^^neiW) , (5) 

where 

n^AM = I d's.Tlis,) Tiib- s,) = I d's.Tlib- S2) Ti{s,) , (6) 

where now T^^ and Ta- are nucleus thicknesses of protons and neutrons, respectively. We 
have assumed one universal inelastic cross section cr'^n^iW) = cr^ei^)- When properly 
normalized, Eq.® describes the probability of the binary af3 collisions and is a simple 
generalization of Eq.(jll). Of course 

TMM{b) = j:n^AM- (7) 

As a second limiting case we shall consider the wounded nucleon model We shall 
assume that the production of particles is proportional to the number of wounded nucleons 
rather than to the number of binary collisions. The simplest version of the model must 
be generalized to include different distributions of protons and neutrons and differences of 
elementary cross sections. Somewhat schematically the cross section for the nuclear collision 
can be written as 



^ a,p=p,n 



(8) 



In the formula above N^°^_ is the number of wounded a (p or n) in nucleus Ai or A2 and 
lyf is the probability that the wounded a interacted with /5 (p or n) from nucleus A2 or Ai, 
respectively. Eq.® is equivalent to Eq.(|T]) with 

N^p-..^ = KJl {b, W) wl{b) + ATJ/X (6, W) w^b) . (9) 

By construction, in our approach the numbers of wounded protons and neutrons repro- 
duce the well known formula from jl^J for the number of wounded nucleons 

N^/l ib) = N-;i (b) + KJl (6) . (10) 



Numbers of wounded protons and neutrons were calculated in the zero-range nucleon-nucleon 
approximation. In calculating the number of wounded nucleons we have taken = 30 
mb, relevant for CERN SPS energies. Our construction requires also 

wf{b)+wnb) = l. (11) 

In this communication the fractions wf were estimated as 

ni'^ih] = ' ' (12) 

which by construction fulfils (lllj). With this prescription Eq.Q reduces to: 

iV„^_± = (iV-/X(6, W) + iV-- (6, W)) w^{b) w^,{b) . (13) 

In the present note we use proton and neutron densities calculated in the Hartee-Fock- 
Bogoliubov method |lo| with Skyrme interaction SLy4 (see Fig^. 

Above we have only sketched the model. More details will be presented and discussed in 
detail elsewhere Q]. 

In Fig 121 we present the ratio -R+/_ (see Eq.Q) as a function of Xp for different values of 
the impact parameter b in the two models considered. For comparison we show preliminary 
experimental data for "central collisions" (solid circles) and "peripheral collisions" (open 
circles) presented in The notion of central and peripheral collisions was not specified 
in Therefore the data can be used only as an indication of the effect. At present no 
quantitative comparison of model results to experimental data is possible. Our approach 
explains the experimental data provided they are extremely peripheral (this requires their 
more detailed analysis). Our results are reliable up to xp ~ 0.2. It is known that the 
HIJING code does not describe the tt^-tt^ asymmetry in elementary collisions in the region 
of large xp. 

For completeness in Figj^we present the ratio -R+/_ as a function of the impact parameter 
b for different values of xp- Surprisingly the binary collision picture gives very similar results 
to the predictions of the wounded nucleon model. 

The predictions of the two models are almost identical. This suggests that a detailed 
comparison of model results with the well defined [xp, b) experimental data could open a 
new possibility to study the neutron density profile. We expect that the NA49 collaboration 
at CERN will be able to gather the corresponding experimental data in the near future. 



Can it provide a method competitive to that offered by proton-nucleus elastic scattering, 
antiprotonic atoms or parity violating electron scattering? Of course results of these methods 
must finally converge. Therefore one may hope that together they will provide more reliable 
information on neutron distribution in nuclei. 
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FIG. 1: Neutron-to-proton ratio jr for "^^^Pb. The experimental point is taken from [J] while 
the experimental grey band from J^. The result of the HFB calculation |1Q] is shown by the thin 
solid line. 
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FIG. 3: The ratio as a function of the Feynman xp for the model with binary colhsion (left 

panel) and for the model with the scaling with number of wounded nucleons (right panel). Also 
shown are the preliminary experimental data from Ref . [l^ . 
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FIG. 4: The ratio as a function of impact parameter for the model with binary colHsion 

(left panel) and for the model with the scaling with number of wounded nucleons (right panel). 
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